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A  guidance  and  control  (G&C)  design  tool  to  rapidly  assess  the  necessary  control  effort 
of  a  conceptual  space  access  vehicle  to  track  its  flight  trajectory  is  described.  This  tool 
can  be  used  as  part  of  the  preliminary  design  cycle  in  configuration,  trajectory  planning, 
structural  analysis,  aerodynamic  modeling,  control  surface  sizing.  Given  a  conceptual  con¬ 
figuration  for  a  space  access  vehicle  and  a  desired  trajectory  for  a  reentry  flight,  this  G&C 
tool  provides  an  inner  loop  feedback  control  law  and  outer  loop  feedback  guidance  law  to 
track  the  given  trajectory.  The  inner  loop  control  law,  based  on  dynamic  inversion  with  a 
non-linear  control  allocator,  is  used  to  linearize  the  vehicle  dynamics  over  its  flight  envelope 
and  assign  control  tasks  over  the  available  control  effectors  to  track  the  desired  roll  rate 
P,  pitch  rate  Q,  and  yaw  rate  R .  The  outer-loop  guidance  law  is  based  on  backstepping 
method  that  transforms  the  trajectory-related  flight  path  angle  7  and  desired  bank  angle  0 
into  commands  in  roll  rate  Pcomd,  pitch  rate  Qcomcz,  and  yaw  rate  RCOmd •  Assessment  of  the 
vehicle’s  tracking  performance  and  associated  aero-control  usage  can  be  made.  This  assess¬ 
ment  can  then  be  used  to  determine  the  appropriate  control  sizing,  and  mass  properties 
for  the  vehicle. 
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I.  INTRODUCTION 


During  the  space  access  vehicle  preliminary  design  process,  it  is  necessary  to  quickly  and  economically 
assess  the  vehicle’s  ability  to  track  a  given  trajectory  utilizing  aerodynamic  and  propulsion  control  effectors. 
This  paper  details  a  framework  that  efficiently  incorporates  the  space  access  vehicle’s  stability,  guidance  and 
control  considerations  into  the  initial  configuration  development.  In  this  approach,  a  well-known,  high-fidelity 
trajectory  generator  (Program  to  Optimize  Simulated  Trajectories),  a  fast  aerodynamic  data  computation 
algorithm  (Missile  Datcom),  and  a  robust,  large  flight  envelope  control  law  are  integrated  in  the  analysis 
and  assessment  process  to  evaluate  vehicle  performance  in  stability,  guidance,  and  control.  POST  (Program 
to  Optimize  Simulated  Trajectories)  is  a  trajectory  computation  program  developed  by  NASA-Langley  in 
1970 ’s  to  support  the  Space  Shuttle  program.  POST  finds  a  user-defined  optimal  trajectory  based  upon 
a  simulation  model  with  performance  and  loading  constraints.  This  optimal  trajectory  is  a  compromise 
between  the  ascent  phase  trying  to  maximize  the  payload  delivered  to  orbit,  the  entry  phase  trying  to  limit 
re-entry  aerodynamic  heating  and  structural  loads,  and  the  approach  phase  trying  to  prepare  the  vehicle 
for  a  successful  landing.  The  ascent  phase  ends  at  stage  separation  which  is  defined  as  a  velocity  of  7,  000 
feet  per  second,  a  flight  path  angle  of  20  degrees,  an  angle  of  attack  of  0  degrees,  and  an  altitude  higher 
than  160,000  feet.  During  this  phase  the  vehicle  is  ascending  rapidly,  passing  through  the  transonic  sound 
barrier,  and  trying  to  reach  stage  separation  criteria.  Flight  loads  are  constrained  or  limited  to  a  dynamic 
pressure  of  less  than  800  pounds  per  square  foot,  body  axial  loads  less  than  6  g’s,  and  body  normal  loads 
less  than  3.5  g’s.  The  entry  phase  begins  at  stage  separation  and  ends  near  the  landing  site  which  is  defined 
to  be  an  altitude  greater  than  15,  000  feet  and  a  velocity  greater  than  Mach  0.3.  During  this  phase  the 
vehicle  typically  experiences  the  peak  heating,  dynamic  pressure,  and  normal  loads  because  the  vehicle  is 
ballistic  and  aerodynamic  controls  are  ineffective.  In  this  phase  the  vehicle  will  descend  at  a  constant  angle 
of  attack  of  35  degrees  and  then  perform  a  pull-up  maneuver  to  intercept  a  constant  flight-path  angle  glide 
of  12  degrees.  The  imposed  constraints  are  the  same  dynamic  pressure,  axial  g,  and  normal  g  loads  cited 
above.  Heating  constraints  are  not  enforced  because  temperatures  estimates  are  calculated  using  MINIVER 
off-line.  The  heat  flux  or  heat  rate  estimates  calculated  in  POST  are  used  for  reference  only  and  ’’ball-park” 
approximation. 

The  POST  simulation  model  consists  of  a  rocket  engine  model,  an  aerodynamic  model,  mass  model, 
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atmosphere  model,  and  Earth  gravity  model.  All  of  these  models  are  combined  to  estimate  the  forces  and 
moments  exerted  upon  the  vehicle  over  time.  Typically,  POST  trajectories  use  three  degrees  of  freedom  to 
solve  the  optimization  problem;  therefore,  losses  due  to  trim  effects  are  not  considered  and  control  effectors 
are  not  sized  appropriately.  In  some  cases  vehicles  may  appear  to  ”  close”  (correctly  sized  in  terms  of  required 
fuel  and  subsystems)  that  are  actually  not  closed  and  may  not  be  closeable  (fuel  and  subsystem  growth  rate 
is  too  high).  Therefore,  accurate  sizing  makes  six  degree  of  freedom  (6  DOF)  simulations  a  vital  part  of 
vehicle  design  process.  Additionally,  the  6  DOF  simulations  provide  higher  fidelity  results  that  can  be  used 
to  impact  other  subsystems  such  as  the  thermal  protection  system,  main  rocket  engine  thrust,  and  main 
propellant  tank  volume. 

POST  finds  a  user-defined  optimal  trajectory  based  upon  a  simulation  model  with  performance  and 
loading  constraints.  This  optimal  trajectory  is  a  compromise  between  the  ascent  phase  trying  to  maximize 
the  payload  delivered  to  orbit,  the  entry  phase  trying  to  limit  re-entry  aerodynamic  heating  and  structural 
loads,  and  the  approach  phase  trying  to  prepare  the  vehicle  for  a  successful  landing.  The  ascent  phase 
ends  at  stage  separation  which  is  defined  as  a  velocity  of  7,000  feet  per  second,  a  flight  path  angle  of  20 
degrees,  an  angle  of  attack  of  0  degrees,  and  an  altitude  higher  than  160,000  feet.  During  this  phase  the 
vehicle  is  ascending  rapidly,  passing  through  the  transonic  sound  barrier,  and  trying  to  reach  stage  separation 
criteria.  Flight  loads  are  constrained  or  limited  to  a  dynamic  pressure  of  less  than  800  pounds  per  square 
foot,  body  axial  loads  less  than  6  g’s,  and  body  normal  loads  less  than  3.5  g’s.  The  entry  phase  begins 
at  stage  separation  and  ends  near  the  landing  site  which  is  defined  to  be  an  altitude  greater  than  15,000 
feet  and  a  velocity  greater  than  Mach  0.3.  During  this  phase  the  vehicle  typically  experiences  the  peak 
heating,  dynamic  pressure,  and  normal  loads  because  the  vehicle  is  ballistic  and  aerodynamic  controls  are 
ineffective.  In  this  phase  the  vehicle  will  descend  at  a  constant  angle  of  attack  of  35  degrees  and  then  perform 
a  pull-up  maneuver  to  intercept  a  constant  flight-path  angle  glide  of  12  degrees.  The  imposed  constraints  are 
the  same  dynamic  pressure,  axial  g,  and  normal  g  loads  cited  above.  Heating  constraints  are  not  enforced 
because  temperatures  estimates  are  calculated  using  MINIVER  off-line.  The  heat  flux  or  heat  rate  estimates 
calculated  in  POST  are  used  for  reference  only  and  ”  ball-park”  approximation. 

The  POST  simulation  model  consists  of  a  rocket  engine  model,  an  aerodynamic  model,  mass  model, 
atmosphere  model,  and  Earth  gravity  model.  All  of  these  models  are  combined  to  estimate  the  forces  and 
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moments  exerted  upon  the  vehicle  over  time.  Typically,  POST  trajectories  use  three  degrees  of  freedom  to 
solve  the  optimization  problem;  therefore,  losses  due  to  trim  effects  are  lost  and  control  effectors  are  not 
sized  appropriately.  In  some  cases  vehicles  may  appear  to  ”  close”  (correctly  sized  in  terms  of  required  fuel 
and  subsystems)  that  are  actually  not  closed  and  may  not  be  closeable  (fuel  and  subsystem  growth  rate 
is  too  high).  Therefore,  accurate  sizing  makes  six  degree  of  freedom  (6  DOF)  simulations  a  vital  part  of 
vehicle  design  process.  Additionally,  the  6  DOF  simulations  provide  higher  fidelity  results  that  can  be  used 
to  impact  other  subsystems  such  as  the  thermal  protection  system,  main  rocket  engine  thrust,  and  main 
propellant  tank  volume. 

The  aerodynamic  data  of  the  vehicle  can  be  obtained  from  Missile  Datcom  which  is  a  tool  to  rapidly 
estimate  the  aerodynamics  of  a  wide  variety  of  vehicle  configurations.  The  predictive  accuracy  of  Datcom  is 
adequate  for  preliminary  designs.  Iterations  on  the  vehicle  configurations  are  inevitable  since  the  ultimate 
shape  of  the  vehicle  will  be  dependent  upon  the  subsystem  being  utilized,  such  as  payload  size,  propulsion 
method,  launch  and  landing  configurations.  Once  the  optimal  trajectory  has  been  calculated  by  POST 
for  a  given  vehicle  configuration  having  the  aerodynamic  properties  indicated  by  DATCOM,  the  task  of 
stabilizing  the  vehicle  and  tracking  the  optimal  trajectory  is  carried  out  by  a  control  law.  The  control  law  is 
designed  for  the  flight  envelope  and  can  accommodate  the  drastic  change  in  speed,  altitude,  and  vehicle’s  mass 
properties  during  ascent  and  reentry.  These  changes  in  the  flight  environment  result  in  substantial  variations 
in  aero-dynamic  pressure,  stagnation  temperature,  center  of  gravity,  and  moments  of  inertia.  Furthermore, 
the  launch  vehicle  is  powered  by  a  propulsion  system  during  its  ascent  and,  for  a  reusable  vehicle,  may  be 
unpowered  during  reentry.  The  vehicle’s  launch  configuration  and  landing  configuration  may  also  be  different 
from  each  other.  The  novel  control  law  seeks  to  stabilize  the  vehicle  and  track  an  optimal  trajectory  without 
the  lengthy  design  process  or  a  complicated  control  law  gain  scheduling  that  is  traditionally  required.  The 
inputs  to  the  inner-loop  control  law  are  the  commanded  roll  rate  PCOmd ,  commanded  pitch  rate  QCOmd ,  and 
commanded  yaw  rate  RCOmd •  The  guidance  and  control  interface  translates  the  bank  angle  command  and 
angle  of  attack  command  into  the  commanded  roll,  pitch  and  yaw  rates  PCOmdi  Qcomdi  Pcomd  respectively. 
Using  the  output  of  the  control  law,  the  vehicle  designer  can  assess  vehicle  performance  in  tracking  the 
desired  trajectory  and  make  modifications  to  its  configuration,  if  necessary.  The  design  process  is  shown  in 
Figure  1. 
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Vehicle  Configuration 
Designer 


Figure  1.  Space  Access  Vehicle’s  Preliminary  Design  Cycle 

II.  Vehicle  Modeling 

The  preliminary  configuration  for  baseline  space  access  vehicle  is  shown  in  Figure  2.  The  body  as  modeled 
by  Missile  Datcom  consists  of  a  blunted  nose  followed  by  a  cylindrical  afterbody  with  a  17  ft  diameter.  A 
single  body  flap  was  modelled  at  the  base  of  the  body.  There  was  no  modelling  of  external  rocket  nozzles. 
A  straked  wing  with  an  outer  panel  sweep  of  42°  was  modelled.  Vertical  tails  were  placed  on  each  wingtip. 
Five  control  devices  were  modelled,  two  rudders  (one  on  each  vertical  tail),  two  elevons  (one  on  each  wing) 
and  a  body  flap.  The  chord  of  the  rudders  was  assumed  to  be  0.30  of  the  local  chord.  The  chord  of  the 
elevons  varied  linearly  from  0.20  of  the  wing  root  to  0.50  of  the  wing  tip.  As  defined,  the  wing  root  is  20.4 
ft  and  does  not  include  the  strake.  Missile  Datcom1  was  used  to  calculate  the  aerodynamic  characteristics 
of  the  vehicle.  Missile  Datcom  is  a  widely  used  engineering  level  code  that  uses  the  component  buildup 
technique  to  predict  vehicle  characteristics.  Code  input  consists  of  body,  wing  and  tail  geometry,  Mach 
number,  altitude,  angle  of  attack  and  control  deflections.  Control  devices  are  limited  to  all  moving  surfaces 
or  plain  trailing  edge  flaps.  At  each  flight  condition  the  six-body  axis  force  and  moment  coefficients  are 
provided.  Both  theoretical  and  empirical  methods  are  included  that  encompass  the  entire  speed  regime  from 
subsonic  to  hypersonic.  Missile  Datcom  has  been  shown  to  provide  very  good  agreement  with  experimental 
data  for  a  variety  of  configurations.  To  validate  the  code  for  RLV  type  configurations,  extensive  comparisons 
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Figure  2.  Preliminary  Configuration  for  A  Space  Access  Vehicle 

have  been  made  with  wind  tunnel  data  for  the  X-34  and  X-40  configurations.  Some  of  the  X-34  comparisons 
are  given  by  Ngo  and  Blake.2 

Moments  of  inertia  for  the  ARES  vehicle  were  calculated  using  the  following  equations  from  Roskam:3 


Ixx  — 

—  (kx  b)2 

g 

(i) 

II 

1? 

to 

(2) 

Izz  — 

w  b  +  l2 
g[kz  2  } 

(3) 

The  non-dimensional  radii  of  gyration  (k  factors)  were  taken  from  an  AFRL  database  of  re-entry  vehicle 
designs  such  as  the  Space  Shuttle,  X-40,  etc..  Values  used  for  kx ,  ky  and  kz  were  0.150,  0.25  and  0.30 
respectively.  This  method  does  not  give  the  product  of  inertia  (Ixz)  so  this  was  assumed  to  be  zero.  The 
dimensions  of  the  baseline  space  access  vehicle  design  is  summarized  in  Table  1.  The  control  effectors  limits 
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are  shown  in  Table  2 


Table  1.  Space  Access  Vehicle  Configuration  Properties 


Fuselage  Length 

Wing  Span 

Weight 

I XX 

lyy 

Izz 

Ixz 

(feet) 

(feet) 

(pounds) 

(slug- feet2) 

(slug-feet2) 

(slug- feet2) 

(slug- feet2) 

102 

42 

65,474 

81,000 

1,322,00 

950,000 

0 

Table  2.  Space  Access  Vehicle  Control  Effectors 


Left  Elevon 

Right  Elevon 

Left  Rudder 

Right  Rudder 

Body  Flap 

=b  30° 

±30° 

±30° 

=b  30° 

±  20° 

III.  DYNAMIC  INVERSION 

The  inner- loop  control  architecture  developed  in  this  work  consists  of  three  major  components:  a  dynamic 
inversion  control  law,  a  control  allocation  algorithm,  and  precompensation.  The  purpose  of  the  inner-loop 
control  system  is  to  accurately  track  body- frame  angular  velocity  vector  commands.  The  goal  of  dynamic 
inversion  in  flight  control  is  to  cancel  the  wing-body-propulsion  forces  and  moments  with  control  effector 
forces  and  moments  such  that  the  vehicle  can  accurately  track  some  desired  commands.  Dynamic  inversion 
control  laws4  require  the  use  of  a  control  mixer  or  control  effector  allocation  algorithm  when  the  number  of 
control  effectors  exceeds  the  number  of  controlled  variables,  or  when  actuator  rate  and  position  limits  must 
be  taken  into  account.  It  is  quite  common  that  the  desired  control  variable  rate  commands  can  be  achieved 
in  many  different  ways  and  so  control  allocation  algorithms  are  used  to  provide  unique  solutions  to  such 
problems.7,5  The  control  allocation  algorithm  is  significantly  improved  by  including  an  intercept  term.6  To 
complete  the  inner-loop,  precompensation  blocks  are  designed  to  produce  the  desired  closed-loop  dynamics. 
For  the  purpose  of  demonstration,  we  develop  a  dynamic  inversion  control  law  for  a  vehicle  with  five  control 
surfaces.  The  control  surfaces  include  two  rudders,  two  elevons,  and  a  bodyflap.  An  outer-loop  control 
system  generates  body- frame  angular  velocity  commands  (pdes,  Qdes ,  r des ),  that  the  inner- loop  dynamic 
inversion  control  system  attempts  to  track.  The  dynamics  of  the  body-frame  angular  velocity  vector  for  this 

7  of  21 


American  Institute  of  Aeronautics  and  Astronautics 


vehicle  can  be  written  as 


u  =  f(u;,  P )  +  g(P ,  8)  (4) 

where  u  =  [p  q  r]T,  p,  q,  and  r  are  the  rolling,  pitching,  and  yawing  rates,  respectively,  P  denotes  measurable 
or  estimable  quantities  that  influence  the  body-frame  states,  and  8  =  (5i,  £2,  •  •  •  ,  8n)T  is  a  vector  of  control 
surface  deflections.  The  vector  P  contains  variables  such  as  angle  of  attack,  sideslip,  Mach  number,  and  ve¬ 
hicle  mass  properties.  The  term  g(P,  8)  includes  the  control  dependent  accelerations,  while  the  term  f(u,  P) 
describes  accelerations  that  are  due  to  the  base-vehicle’s  (wing-body-propulsion)  aerodynamic  properties. 
The  moment  equations  for  a  vehicle  in  the  body-frame7  can  be  manipulated  to  form  control  dependent  and 
control  independent  terms.  It  is  assumed  that  the  mass  properties  of  the  vehicle  under  consideration  are 
constant,  thus,  the  time  derivative  of  the  inertia  matrix  can  be  set  to  zero,  i.e.,  1  =  0.  Then,  Equation  4  can 
be  written  as 

w  =  r1(GB(w,P,(5)-wxIw)  (5) 


where 


G  B(w,  P,  6)  =  GBAE(u>,  P)  +  Gi(P,  S ) 


- 

- 

L 

L 

M 

+ 

M 

N 

BAE 

N 

(6) 


In  Equations  5  and  6,  I  is  the  inertia  matrix  and  L,  M,  and  N  are  the  rolling,  pitching,  and  yawing  moments. 
In  Equation  6,  G^(w,  P)  is  the  moment  generated  by  the  base  aerodynamic  system  (wing-body-propulsion 
system)  and  G$(P,<5)  is  the  sum  of  moments  produced  by  the  control  effectors.  Therefore, 


(7) 


f  (w,  P)  =  I-1  [GBae(w,  P)  -  w  x  Iw] 
g(P,«5)  =  I-1G5(P,<5) 

In  order  to  utilize  a  linear  control  allocator,  it  is  necessary  that  the  control  dependent  portion  of  the 
model  be  linear  in  the  controls.  Hence,  an  affine  approximation  is  developed  such  that 


G5(P,<5)  «G5(P)<5  +  e(P,d) 


(8) 
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The  term  e(P,  S)  is  an  intercept  term6  for  the  body- axis  angular  accelerations  which  is  used  to  improve  the 
accuracy  of  linear  control  allocation  algorithms.  Using  Equations  4,  7,  and  8,  the  model  used  for  the  design 
of  the  dynamic  inversion  control  law  becomes 

iv  =  f(u),  P)  +  I_1Gi(P)<5  +  I-1e(P,  <5)  (9) 


The  objective  is  to  find  a  control  law,  that  provides  direct  control  over  cb,  so  that  cb  =  cb^es.  Hence,  the 
inverse  control  law  must  satisfy 


wdes-f(u;,P)-I-1e(P)<5)  = 

r'G^p)* 


(10) 


Equation  10  provides  the  dynamic  inversion  control  law  for  the  body-frame  angular  velocity  vector. 


IV.  CONTROL  ALLOCATION 


Since  there  are  more  control  effectors  than  controlled  variables  and  the  control  effectors  are  restricted  by 
position  and  rate  limits,  a  control  allocation  algorithm  is  necessary.  For  the  lifting  body  under  consideration, 
there  are  three  controlled  variables,  namely,  roll,  pitch,  and  yaw  rates,  while  there  are  five  control  surfaces 
(left  and  right  rudders,  left  and  right  elevators,  and  a  bodyflap).  Hence,  a  control  allocation  scheme  must 
be  used  to  insure  that  Equation  10  is  satisfied. 

Control  allocators  are  used  in  conjunction  with  some  type  of  feedback  control  law  whose  output  consists  of 
one  or  more  pseudo-control  commands  (typically  desired  moment  or  acceleration  commands).  The  number  of 
pseudo-control  commands  is  always  less  than  or  equal  to  the  number  of  control  effectors.  Dynamic  inversion 
control  laws  and  control  allocation  algorithms  fit  together  quite  naturally  since  the  pseudo-control  commands 
are  easily  identifiable.  Also,  it  is  quite  common  that  the  desired  commands  can  be  achieved  in  many  different 
ways  and  so  control  allocation  algorithms  are  used  to  provide  unique  solutions  to  such  problems. 

To  begin  development  of  the  allocator,  rewrite  Equation  10  as 

ddes  =  Wries  -  f(w,  P)  -  I_1e(P,  (5)  = 


(11) 


I_1G5(P)<5  =  B<5 

where  d des  are  the  body-axis  accelerations  that  need  to  be  produced  by  the  control  effectors  and  B  is  the 
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control  effectiveness  matrix  defined  as 


8L  8L  m  m  m  8L 
88\  882  88n 


b  =  i-1g5(p)  =  i-1 


8M  8M  m  m  m  8M 
88i  882  88 n 


8N  8N 

88 1  882 


8N 

88n 


The  control  allocation  objective,  in  the  linear  case,  is  to  find  8  such  that 


dries  =  B<5 


(12) 


(13) 


subject  to  rate  and  position  limits  on  the  control  effectors.  Notice  that  Equation  13  defines  a  linear  subspace 
in  the  (d^es,^)  space. 

Equation  13  can  be  posed  as  the  following  optimization  problem: 


min  Je  =  min  ||B5  -  d<fea|li  (14) 

8  8 


subject  to 


6  <5  <5 


(15) 


where  Je  is  the  performance  index  for  the  error  minimization  problem,  8 ,  (5  are  the  most  restrictive  lower 
and  upper  limits  on  the  control  effectors,  respectively  and  the  1-norm  is  selected  so  that  linear  programming 
techniques  can  be  used  to  solve  the  problem.8  More  specifically, 


S  =  min(<5[/,  8  +  8  At) 
5  -  ma x(Sl,  S  —  5 At) 


(16) 


where  8l,Su  are  the  lower  and  upper  position  limits,  <5  is  the  current  location  of  the  control  effectors,  <5  is 
a  vector  of  rate  limits,  and  At  is  the  timestep  or  control  update  rate. 

If  sufficient  control  authority  exists  such  that  Je  can  be  made  identically  equal  to  zero,  then  it  may  be 
possible  to  optimize  a  sub-objective.  This  optimization  problem  can  be  posed  as  follows: 


min  Jc  =  min  ||W5(5  -  cSp)^ 
8  8 


(17) 
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subject  to 


—  d  des 


8  =  minify,  8  +  8  At)  (^) 

8  =  ma x(8l,  8  —  8 At) 

where  WA  is  a  weighting  matrix  and  8P  is  a  preferred  set  of  control  effector  deflections.  The  problem  posed 
in  Equation  17  is  termed  the  control  minimization  problem. 

In  practice,  the  two  optimization  problems  given  in  Equations  14  and  17  are  combined  to  form  what  is 
known  as  the  mixed  optimization  problem.  The  mixed  optimization  problem  is  defined  as 

min  JM  = 

5  (19) 


min  (||B<5  -  +  A  ||W4(5  -  <*p)  Hj 

0 

where  the  parameter  A  is  used  to  weight  the  error  and  control  minimization  problems.  For  this  work,  it  was 
determined  that  A  =  0.01  provided  good  error  minimization  while  still  driving  the  control  effectors  to  the 
preferred  values  when  sufficient  control  authority  existed.  The  advantage  of  the  mixed  optimization  problem 
is  that  it  can  often  be  solved  faster  and  with  better  numerical  properties  as  compared  to  sequentially  solving 
the  error  and  control  minimization  problems.5 


A.  Control  Allocation  Preference  Vector  and  Effector  Failures 

As  specified  in  Equation  19,  a  preference  vector,  8P,  must  be  selected.  One  difficulty  with  the  linear  pro¬ 
gramming  framework  for  solving  the  control  allocation  problem  is  that  no  model  of  the  control  allocator 
exists.  This  causes  problems  when  performing  linear  stability  analysis  as  there  is  no  way  to  know  the  in¬ 
put/output  relationship  of  the  allocation  algorithm.  Fortunately,  when  sufficient  control  authority  exists, 
the  allocation  algorithm  will  attempt  to  minimize  the  difference  between  the  control  deflections  and  a  pre¬ 
ferred  set  of  control  deflections.  One  obvious  choice  for  preference  vector  is  the  pseudo- inverse  solution.  In 
this  case,  when  sufficient  control  authority  exists,  the  control  allocation  algorithm  will  drive  the  surfaces 
to  the  pseudo-inverse  solution.  Hence,  in  a  robustness  analysis,  the  control  allocator  can  be  replaced  by 
the  pseudo-inverse  solution  (assuming  sufficient  control  authority  exists).  The  pseudo-inverse  solution  is  the 
two-norm  solution  to  the  control  allocation  problem  and  can  be  formulated  as  follows: 

min  -  (8  +  c)T W (8  +  c)  (20) 

8  2 
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subject  to 


B5  =  d 


des 


(21) 


where  W  is  a  weighting  matrix  and  c  is  an  offset  vector.  To  solve  this  problem,  first  find  the  Hamiltonian 
(H)  such  that 


H  =  i(5rW(5  +  bTW<5+ 

\sT Wc  +  Vwc  +  £(B S  -  ddes ) 


(22) 


where  £  G  Mn  is  an  as  yet  undetermined  Lagrange  multiplier.  Taking  the  partial  derivatives  of  H  with  respect 
to  S  and  £,  setting  these  expressions  equal  to  zero,  and  rearranging,  gives 


§-W*+i(c*wf- 


1 


-Wc  +  (£B)J  =0 


TVT 


W  6  =  — Wc  —  B  £ 


and 


OH 

~d£ 


=  B£  -  d, 


des 


o 


=>  B<5  =  ddes 
=►  BW-‘W  S  =  ddes 

Substituting  Equation  23  into  Equation  24  yields 


(23) 


(24) 


BW‘[-Wc-B'f]  =  d 


des 


(25) 


Solving  for  £T  in  Equation  25  yields 


e  =-(BW-1BJ)-1[ddes  +  Bc 


(26) 


Substituting  Equation  26  into  Equation  23  produces 


W  6  =  -Wc  +  BT(BW“1BT)-1[ddes  +  Be 


(27) 


Simplifying  Equation  27  gives  the  desired  result 


6  =  Sp  = 


-c  +  W-1BJ  (BW_1BJ  j-1  [ddes  +  Be] 


(28) 
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Equation  28  gives  the  pseudo- inverse  solution.  It  should  be  noted  that  if  an  effector  is  offset,  two  items  must 
be  taken  into  account,  position  offset  (— c)  and  the  moments  generated  by  the  offset  (Be).  For  the  specific 
usage  of  the  pseudo-inverse  control  allocation  solution,  the  weighting  matrix  was  selected  to  be  diagonal, 
such  that, 


W  =  diag 


(29) 


[Wi™  WSlf  WSrr  WSlr  WSsB  W5bf] 


where  ’diag’  represents  a  diagonal  matrix  with  the  entries  along  the  main  diagonal  being  the  weights  asso¬ 
ciated  with  each  control  effector. 

This  control  allocation  formulation  allows  one  to  simulate  a  control  effector  failure  rather  easily.  A  failure 
is  introduced  by  simply  setting  the  lower  and  upper  positions  limits  on  the  effected  control  surface  equal  to 
each  other.  For  a  failed  control  surface,  its  effects  must  also  be  accounted  for  in  the  pseudo-inverse  preference 
vector,  which  requires  two  modifications.  First,  the  location  of  the  failure  must  be  inserted  into  the  offset 
vector.  Here,  the  appropriate  component  of  c  is  set  to  the  negative  of  the  failure  position.  Second,  the 
appropriate  entry  in  the  weighting  matrix,  W,  needs  to  be  increased.  Nominally,  the  entries  in  W  are  one 
and  an  increase  in  the  value  will  place  more  penalty  on  usage  of  that  particular  surface. 


V.  Guidance  Loops 

The  guidance  loops  developed  in  this  work  have  as  inputs  angle  of  attack  error  ae  and  euler  angle  phi 
error  </>e  and  provide  as  outputs  the  commanded  roll,  pitch,  and  yaw  rates.  A  backstepping  approach  is  used 
to  move  from  the  angle  of  attack  and  phi  loops  to  the  body- axis  rate  loops. 

To  derive  the  pitch  rate  command,  the  governing  equation  of  motion  is 

7  =  0 -a  (30) 

Therefore, 

a  =  9  —  7  (31) 

By  definition, 

6  =  q  cos  </>  —  r  sin  </>  (32) 
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and 


L  9  ,~o\ 

7  =  S^“fcost  (33) 

where  L  is  the  total  vehicle  lift,  m  is  the  mass,  V  is  the  velocity,  and  g  is  the  acceleration  due  to  gravity. 
Substituting  Eqs.  32  and  33  into  Eq.  31  yields 


,  -A  L  ,  9 
a  =  q  cos  (p  —  r  sm  cp - —  +  77  cos  7 

mV  V 


(34) 


Prom  dynamic  inversion,  the  pitch  rate  command  becomes 


qdes  =  sec  0  (  ades  +  r  sin  <f>  +  -E  -  ^  cos  7 


(35) 


The  desired  angle  of  attack  dynamics  are  defined  by  the  following  proportional- integral  control  on  a  error: 

kj 


&des  —  (  kpa  T  I  (^des  &) 


(36) 


For  the  lateral  channels,  the  roll  rate  command  is  simply  the  result  of  a  proportional-integral  operation 
on  phi  error.  Thus, 


Pdes  —  (  kpfj,  +  )  i^Pdes  0) 


(37) 


From  the  coordinated  turn  equations,  the  yaw  rate  command  is  computed  using 


rdes  =  ptana  +  —  sine 


(38) 


Equations  37,  35,  and  38  define  the  commanded  or  desired  body  axis  rates  generated  from  errors  between 
the  actual  and  desired  angle  of  attack  a^es  and  desired  roll  angle  <fides- 


VI.  RESULTS 

In  this  section,  we  will  use  the  above  aforementioned  methods  to  obtain  a  guidance  and  control  law  for 
a  space  vehicle  undergoing  a  preliminary  design  evaluation.  Given  a  reentry  trajectory,  a  POST  simulation 
outputs  the  desired  angle  of  attack  ades  and  desired  roll  angle  (j)des  as  part  of  its  calculations.  A  guidance 
law  based  on  the  backstepping  method  discussed  in  section  V  is  designed  to  convert  the  desired  angle  of 
attack  ades  and  desired  roll  angle  (pdes  into  the  commanded  roll  rate  PCOmd ,  the  commanded  pitch  rate 
Qcomd  and  the  commanded  yaw  rate  RCOmd •  The  commanded  body- axis  rates  PCOmd ,  Qcomd >  Rcomd ,  h1  turn, 
are  converted  into  required  control  deflections  by  a  dynamic- inversion  control  law  mentioned  in  section  III. 

14  of  21 


American  Institute  of  Aeronautics  and  Astronautics 


A  control  allocator  based  on  the  discussion  in  section  IV  assigns  the  required  control  deflections  over  the 
control  effectors  according  to  their  available  capabilities.  The  tracking  performance  of  the  control  law  is 
shown  in  Figure  3.  Tracking  errors  PCOmd  ~  Factual ,  Qcomd  ~  Q actual,  Rcomd  ~  Factual  are  small.  These 
small  errors  are  desired  and  expected  since  our  design  method  is  formulated  to  directly  track  PCOmd ,  Qcomd, 
RCom •  The  simulation  performed  in  POST  does  not  include  the  impact  of  additional  drag  associated  the 
control  effectors  activities  in  trimming  the  vehicle  and  tracking  the  given  trajectory,  differences  between 
the  actual  roll,  pitch,  and  yaw  angles  and  those  calculated  by  POST  as  shown  in  Figure  4.  In  Figures  5 
and  6  the  deflection  activities  of  the  left  elevon,  right  elevon,  left  rudder  and  right  rudder  are  reasonable 
without  exceeding  its  rate  and  position  limits.  It  is  noted,  however,  that  the  pitch  flap  briefly  saturates  in 
its  deflection  limit.  For  a  given  trajectory,  it  is  not  desirable  for  any  control  effector  to  experience  either 
a  rate  or  a  position  saturation  under  a  nominal  condition.  To  avoid  control  saturation,  there  are  many 
options  ranging  from  hardware  considerations,  to  software  modifications,  to  trajectory  planning.  In  the  next 
discussion,  we  will  examine  how  the  issue  can  be  addressed  from  the  point  of  view  of  stability  and  control 
namely  control  sizing,  and  center  of  gravity  location. 

A.  Aero-Control  Surface  Resizing 

As  seen  in  Figure  5  the  body  flap  experiences  saturation.  By  enlarging  this  control  effector,  it  is  hoped  the 
saturation  can  be  avoided.  For  this  study,  the  cord  of  the  body  flap  is  increased  by  40%  from  its  baseline 
of  5  feet  to  7  feet.  The  control  activities  of  the  body  flap  is  seen  in  Figure  8  The  body  flap  still  experiences 
saturation.  At  this  stage  of  the  design,  other  changes  on  the  vehicle  can  be  experimented  to  affect  the 
vehicle’s  stability  and  control  such  as  moving  the  vehicle’s  center  of  gravity,  reducing  its  strakes. 

VII.  CONCLUSIONS 

In  this  work,  a  rapid  assessment  tool  for  space  access  vehicle  configurations  in  guidance  and  control 
performance  was  presented.  The  reentry  trajectory  was  found  using  a  well-known,  high-fidelity  trajectory 
generator.  To  track  this  trajectory,  an  inner-loop  control  law  was  designed  for  a  reentry  vehicle  with  five 
control  surfaces.  The  control  law  utilized  a  dynamic  inversion  controller  and  a  linear  programming  based 
control  allocation  algorithm.  Evaluation  of  the  vehicle’s  guidance  and  control  performance  showed  that  the 
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P,Q,R 


Figure  3.  Pitch  Rate,  Roll  Rate,  and  Yaw  Rate  Tracking 

body  flap  experienced  deflection  saturation  during  a  nominal  reentry  flight.  It  is  thus  necessary  to  consider 
modifications  of  the  vehicle  such  as  control  surface  sizing,  moving  the  vehicle’s  center  of  gravity,  resizing 
the  wing,  strake  areas.  Using  the  tool  detailed  in  this  paper,  such  modifications  of  the  vehicle  can  be 
accommodated  during  the  early  design  stages  to  incorporate  the  guidance  and  control  requirements 
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Figure  5.  Left  Elevon,  Righ  Elevon,  Body  Flap  Deflections 
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Figure  6.  Left  Rudder,  Righ  Rudder  Deflections 
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Figure  7.  Velocity,  Angle  of  Attack,  Side  Slip  Angle,  and  Flight  Path  Angle  Tracking 
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Figure  8.  Left  Elevon,  Righ  Elevon,  Body  Flap  Deflections 
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